Clock-transport synchronisation for neutrino time-of-flight measurements 
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I A method to synchronise, at the sub-nanosecond level, clocks used for neutrino time-of-flight 

■ measurements is proposed. Clocks situated near the neutrino source and target are compared 

with a moveable clock that is transported between them. The general-relativistic theory of the 
I procedure was tested and verified in an experiment performed by Hafele and Keating in 1972. 

11^ It is suggested that use of such a synchronisation method may contribute to a precise test of 

the Sagnac effect — a measured velocity greater than c — for neutrinos of the proposed LBNE 
r~| ! beam between Fermilab and the Homestake mine. 
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Q I The aims of the present paper are twofold, firstly to present an analysis of clock transport 

O I synchronisation as a possible calibration method for time-of-fiight measurements in two 

existing neutrino beams [1, 2] and one proposed one [3] and secondly to point out the 
interest of such measurements, at the level of Ins accuracy, as a test of the Sagnac 
CN ■ effect for neutrinos — that the measured velocity is not expected to be exactly equal 

to the speed, c, of light in vacuum [4, 5]. The clock transport synchronisation method 
is not proposed as a stand-alone one but rather as a cross-check of continous methods 
^ i such as GPS Common- View GPS [6] as used for the published neutrino time-of-fiight 

d I measurements or Two Way Satellite Time Transfer [7] which have a comparable precision 

Ins or better) than the method proposed here. 

In the region of the Earth, the proper time interval, dr, of a clock is related to the 
interval of 'coordinate time', dt, by the Schwartzschild metric equation [8, 9]: 



dr 



20E 1 / ,2,2 



dt (1) 



where Vr, ve and are components of the velocity of the clock in a spherical polar co- 
ordinate system with origin at the centre of the Earth, and polar axis in the south-north 
direction, fixed in the Earth Centered Inertial (ECl) frame. The ECl frame is a non- 
rotating inertial frame instantaneously comoving with the centroid of the Earth. The 
cordinate axes of the ECl frame therefore have fixed directions relative to the Celestial 
sphere. The quantity 0e is the gravitational potential of the Earth, which, on the as- 
sumption of a spherical Earth is given as: 0e = —GMe/t at distance r from the centre 
of the Earth, where Me is the mass of the Earth. As can be seen by setting {7=0, 
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0E = in (1), coordinate time is the proper time of any clock at rest in the ECI frame 
sufficiently far from the Earth that all gravitational effects of the latter may be neglected. 
In applying (1) it is assumed that the spatial variation of the gravitational fields of the 
Sun, Moon and other members of the Solar System, that give rise to tidal effects, may be 
neglected, as well as the rotation of the Earth around the Sun that changes the angular 
velocity of a rotating. Earth-fixed, frame relative to the ECI frame. These are the same 
approximations that are made in apphcations of (1) to clocks in the sateUites of the GPS 



system [10, 11]. 

For a clock at rest on the Surface of the Earth Vr — — so that (1) simplifies to: 



where R^y is the distance of the clock from the axis of rotation of the Earth and is the 
angular velocity of the latter relative to the ECI frame. Due to its rotation, the Earth 
has the form of an oblate spheroid, the surface of which is known as the geoid. Locally, 
in the oceans, the plane of the geoid is parallel to the surface of the sea. 

It was shown by Cocke [12] that the distortion of the Earth from a spherical shape 
modifies the external gravitational potential of the Earth in such a way that, to a very 
good approximation, the quantity 20e — ^^Rly in (2) is constant for all points on the 
geoid so that identical clocks at the poles or on the Equator will be observed, from the 
ECI frame, to run at the same rate". 

If a clock moves, at fixed r, relative to the surface of the Earth, with speed v, the 
elapsed proper time interval At is given by integrating (1): 



where is the speed, in the ECI frame, of an object, at the same position as the clock, 
but at rest relative to the Earth, and where, in the last member, only terms linear in 
0e/c^ and quadratic in v/c are retained. 

As suggested by Hafele [15, 16] the correctness of the relation (3) was verified in an 
experiment performed by Hafele and Keating (HK) [17] in 1972. An array of caesium- 
beam atomic clocks was flown around the Earth in commercial airliners in West to East 
(W-E) and East to West (E-W) directions. They were compared, before and after the 
flights, with reference Earth-bound clocks at the U.S. Naval Observatory. In agreement 
with prediction, time interval differences, relative to the reference clocks, of 273ns for the 
W-E and -59ns for the E-W flights were observed [17]. 

In the apphcation of Eq. (3) suggested, in the present work, to synchronise clocks 
for time-of-flight measurements, comparisons are made between the epochs recorded by 
similar clocks (i.e. ones running at the same frequency), one at the position of the source, 
and the other at the position of the detector, of the particles. If the clocks are placed 
on the geoid they will be observed from the ECI frame to run at the same rate. As 
will be shown, the clocks can be synchronised by comparing their epochs with that of a 
third similar clock that is transported between them. The single fixed reference clock of 

"Not as stated in Einstein's 1905 special relativity paper [13], slower at the equator. See [14] for a discussion of 
this point. 
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the HK experiment is therefore replaced by two similar clocks at different, but known, 
spatial locations. In order to make this comparison the proper time intervals Atq, At 
recorded by the Earth-fixed and moving clocks, respectively, for the same interval, At, of 
coordinate time are calculated using Eq. (3) taking account of the variation of the values 
of 0E, Vet and v along the path followed by the moveable clock. For the purposes of the 
calculation presented here it is assumed that the Earth is spherical; in actual experimental 
applications the actual shape of the geoid [11] should be taken into account. 

In a practical realisation of the synchronisation method suggested here the clock might 
be equipped with a GPS receiver and associated hardware to record the precise position 
of the clock in the Earth-fixed frame at known epochs during the transport so that the 
integral in Eq. (3) can be accurately evaluated. In the present work the dependence of 
At on the parameters of the clock trajectory is studied within a simple model where 
only the mean speed and altitude of the transport are considered. This is sufficient to 
quantitatively demonstrate the predictions for Ar given by different choices of transport 
trajectory for the three neutrino time-of-ffight experiments that are considered. 

Considering the specific case of neutrino time-of-fiight measurements in the CNGS 
neutrino beam [2, 22] the proper time intervals recorded by reference clocks at CERN 
(At^) and Gran Sasso (At^^), assumed to lie on the geoid, are given by (3), to first 
order in 0e and second order in (3 — vq/ c as 



Ath 



Ato^ 



Ato«s 



/ 



1 + 



1 



dt 



1 + 



dt 



(4) 



since 4>e — {yuY /2 is invariant on the geoid. The proper time interval At recorded 
by a clock transported from CERN (C) to Gran Sasso (GS) during the interval Ai of 
coordinate time is, from (3) and (4): 



At = Ato + 



dt 



(5) 



where h is the altitude above the Earth's surface. Using a mean value approximation: 



f{t)dt {fit)) / dt 



to evaluate the integral in (5), and assuming transport at constant v, gives: 



At - Ato At - Atq _ St GM^{h) 



At 



Ato 



Ato 



'i^^f + ^^ + 2(i;a)-^ (6) 



Since the values of Sr in (6) are of the order of a few nanoseconds precise modelling of 
the altitude and speed of the travelling clock in order to evaluate the integral in (5) is 
not required. For example a 5% uncertainy in the evaluation of the integral gives only 
0.05ns uncertainty in Sr when 6t = 1 ns. 

Introducing a primed. Earth-centered, Earth-fixed, polar coordinate system with north- 
pointing polar axis, (f)' — at the prime meridian and A' the angle of latitude, the coor- 
dinates of CERN are [18]: X'q = 46.05°, 0'c = 6.08° and of Gran Sasso are: A'^g = 42.3°, 
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= 13.58°. If and k' are unit vectors along the x' (A' = 0, 0' = 0), y'{X' = 0, 

0' = 90°) and z' axes of the corresponding right-handed Cartesian coordinate system, the 
spatial separation of the neutrino source at CERN from the OPERA detector at Gran 
Sasso is d = -^(Aa;')^ + (At/')^ + (A^;')^ where: 

Ax' = i?(cos Xqs cos - COS A'c cos 0c) = 0.02936i? (7) 
Ay' = i?(cos A'cg sin - cos A'c sin 0c) = O.lOOSi? (8) 
Az' = i?(sinA'Gg -sinA'c) = -0.0476i? (9) 

The velocity vector v parallel to the neutrino beam is then: 

V = 't;(0.2557£' + 0.8734/ - 0.4144^') = vv (10) 

Also 

(vn) = nR{-i' {cos X' sin (I)') + /(cos A' cos 0')) = Qi?(-0.1236z' + 0.7047/) (11) 

so that 

(v^) ■ V ^ 0.5839nRv (12) 

and 

{(ij^)^) = n^R^icos^x') ^o.Buen^R^ (13) 

Substituting (12)-(14) in (6) with QR = 465 m/s gives: 

5t = 2:221^ ri9.57(/i) - 7157 - - 543d ns CERN - CNGS (15) 

V 

with (/i) in m and v in m/s. Here the relation 

(iCcs 730.53 x 10=^ x 1.0022 732.14 x 10=^ 
At ~ Ato = = = 

V V V 

is used. The straight line distance d from CERN to the OPERA detector at Gran 
Sasso [18] is corrected by the factor Cqs which accounts for clock transport along a 
Grand Circle. 

A similar calculation for the MINOS experiment [1] in the Fermilab (FNAL) NuMI 
beam with d = 732.7 km, Cqs = 1.0022, gives: 

St = fL^^f^ [l9.57(/i) + 11222 - + 2S0v] ns FNAL - MuMI (16) 

Here FNAL is assumed to be at latitude: Ap = 41.85°, longitude: 0f = -88.29° and the 
MINOS detector in the Soudan mine at A^. = 47.81°, 0'g = -92.24°. 

For the proposed long-baseline FNAL-Homestake neutrino beam [3] LBNE, with d — 
1298 km, Cgs = 1-0069 it is found that: 

St = ^-^^^ [l9.57(/i) + 4670 - + 657v] ns FNAL - LBNE (17) 
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Figure 1: Procedure to synchronise clocks A and B at different spatial locations on the geoid by 
comparison with a clock, M, transported at constant velocity between them. See text for discussion. 
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Table 1: Values of vq and Stq at which 6t has a maximum value, as a function of {h) , for 
CERN-CGNS neutrino beam. 
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Table 2: Values of St in ns as a function of v and (h) for the CERN-CNGS neutrino beam. 

where the Homestake mine at Lead, South Dakota, has coordinates: = 44.33°, 0^ = 
-103.83° [19]. 

The procedure to synchronise a clock, ^4, at the source laboratory, registering an 
epoch ta with clock S, at the detector laboratory, registering an epoch tb is as follows. 
Before the clock transport, the difference, Ar^, between the epochs registered by clock 
A and the moveable clock M are recorded at some instant 1. Ar^ = t^^^ — t^^^ . After 
the clock transport the difference between the epochs of B and M is recorded at some 
instant 2: Ar^ = Tb 'm- -'-^ order to synchronise clock B with clock A the quantity 
At = Atb — Ata -I- St (see Fig. 1) must be subtracted from the epoch registered by clock 
B. Note that if all three clocks have the same frequency when at rest in the same frame of 
reference the differences of epochs Ar^ (Ar^) may be recorded at any time before (after) 
the transport of M, so that the instants 1 and 2 are arbitary. 

Curves of 6t versus v for various values of (h) as calculated by (15), (16) and (17) 
are shown in Figs. 2, 3 and 4 respectively. A feature of all three neutrino beams is 
the existence, for some values of (h), of 'magic' values vm of the velocity for which St 
vanishes (given by the intersection of the curves with the dashed line) and so gives the 
possibility of a simple and accurate synchronisation procedure, as if the moving clock 
was equivalent to the ground-based ones. Magic velocity values correspond to exact 
cancellation of the spccial-rclativistic time dilation effect, which slows the rate of the 
moving clock as observed in the ECI frame, and the general-relativistic gravitational blue 
shift which increases its rate. Curves of vu versus (h) for all three beams are shown in 
Fig. 5. For values of v: 50 km/h < v < 1000 km/h corresponding to road, rail or air 
transport St takes, in all cases, absolute values of a few nanoseconds or less. 

The geometry of the CERN-CNGS beam is such that, for small values of (h), the 
curves of St versus v have maximum values. The corresponding velocity vq and time 
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Figure 2: Curves of 5t versus v for different values of (h) for the CERN-CNGS neutrino beam. 
Magic values of v are given by the intersections of the curves with the dashed line St — 0. 
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Figure 3: Curves of St versus v for different values of (h) for the FNAL-MuMI neutrino beam. 




interval Stq are given by the solution of the equation d{ST)/dv — as: 



^^0 



_ ,2, ,^.2 2GMe(/i) 



{{vnr) - {i^r - 



(18) 



5to = -^^{vo + v-{v^)) (19) 

The largest value of {h) for which a maximium value of 5t exists is ~ 365 m corre- 
sponding to vq — Q. The dependence of vq and 5tq on {h) is given in Table 1. Over the 
range m < {h) < 300 m, Stq varies by only ±0.2ns about an average value of 2.7ns. As 
shown in Table 2 the feature of a weak dependence of 6t on both v and (h) persists for 
365 m < (h) < 1000 m where the curves are monotonic. For example, for v =500 km/h 
and {h) = 500m one percent uncertainties in v and {h) give only ~ 0.24 % and ~ 0.09 % 
uncertainties, respectively in 5t, suggesting that this quantity may be determined with 
a precision of one tenth of a nanosecond or better. In order to achieve this precision 
over the typical beam length of ~ 1000 km, with 100 km/h < v < 1000 km/h requires a 
frequency stability of the clock in the range 3 x 10^^^ — 3 x 10^^'^/h, compatible with the 
typical frequency stability of a precision atomic clock of Ins/day. 

Since the Great Circle route from CERN to Gran Sasso passes over the Alps a possible 
strategy to keep h below 1 km to take advantage of the almost constant value of 5t as 
a function of v and (/i), as shown in Table 2, is to transport the clocks by air above the 
Rhone to the Mediterranean, then by a Grand Circle route to Gran Sasso. An alternative 
strategy is a Grand Circle route directly over the Alps at an altitude of 5- 10km at or 
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Figure 5: Magic velocity values as a function of (h) for different neutrino beams. 

near the magic value of v, again enabling a precise synchronisation of the local clocks at 
CERN and Gran Sasso. 

For the FNAL neutrino beams: 



Inspection of Eq. (18) shows that no real value of Vq exists for these beams; the curves of 
St versus v are therefore monotonic for all values of (h). As can be seen in Figs. 3 and 
4 magic values of v are obtainable by conventional air transport at speeds of 1000-2000 
km/h for the MuMI beam but not for LBNE that requires speeds in excess of 2000 km/h. 

It is interesting to compare the uncertainty in the clock synchronisation constant St 
with the variation of neutrino times-of-fiight due to the Sagnac effect [4, 5, 20, 21], i.e. the 
effect of the motion of the neutrino target, due to the Earth's rotation, on the observed 
time-of- flight. For neutrinos with energy much greater than their mass, the relative speed 
Cr of the neutrino and the target is, at lowest order in the velocity if^ of the target, given 



where v is a, unit vector in the direction from the neutrino source to the target. For the 



{{vnf) - {v^f = -11.22 X 10^ MuMI 
= -4.67 X 10^ LBNE 



by [4, 5]: 




(20) 
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three neutrino beams considered above it is found that: 



-9.05 X 10 



.-7 



CNGS 



c 



= 4.65 X 10 



,-7 



MuMI 



= 1.10 X 10 



-6 



LBNE 



The corresponding difTerences of time-of-flight compared with propagation at speed c are: 
2.2 ns [4], -1.14 ns and -4.74 ns respectively, to be compared with the possible uncertainty 
in 8t of the order of a nanosecond or less. Especially the proposed FNAL-LBNE beam [3] 
is well adapted to a precise test of the Sagnac effect for neutrinos. 

In the case of four experiments which have published neutrino time-of-flight measure- 
ments: MINOS [1], OPERA[22], ICARUS[23] and LVD [24] the dominant systematic 
errors were related to local measurements of the times of production and detection of the 
neutrinos rather than errors in the synchronisation of the precision master clocks near to 
the source and target. In the LVD measuremet that has the smallest quoted systematic 
uncertainty of 3.3 ns on the time-of-flight measurement, only 1 ns is assigned to GPS 
synchronisation. In order to be sensitive to the Sagnac effect more accurate local timimg 
of neutrino production and detection events than hitherto achieved is required in addition 
to synchronisation uncertainty at the 1 ns level or better. 
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